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METHOD FOR DETERMINING THE INTERSECTION 
OF POLYGONS USED TO REPRESENT GEOGRAPHIC FEATURES 

BACKGROUND OF THE INVENTION 

The present invention relates to a method of determining the intersection of two 
polygons, particularly with respect to determining the intersection of two polygons that 
represent geographic features. 

In geographic databases, polygons are used to represent two-dimensional areas, 
such as lakes, recreational parks, cities, counties, states, golf courses, and so on. A 
commonly needed function performed using a geographic database is to find that part of 
one two-dimensional geographic feature (such as a lake, a recreational park, etc.) which 
is located inside another two-dimensional geographic feature (such as a city, county, 
state, and so on). The solution can be found by determining the intersection of the 
polygons used to represent the two-dimensional geographic features. 

There are prior methods for determining the intersection of polygons. One well 
known algorithm for determining the intersection of polygons is described in 
Computational Geometry (O'Rourke, Joseph. Second Edition. Cambridge University 
Press, 1998). This algorithm takes all shape points from two polygons, sorts them by 
their y coordinate value and sweep searches the line segments at each shape point from 
top to bottom. If a line segment is located in both polygons, then it is an intersection 
boundary. With this prior algorithm, if there are n shape points in two polygons, the 
overall complexity is n*log(n). Thus, although this prior algorithm works satisfactorily, 
it is computationally intensive and accordingly takes a relatively large amount of time to 
obtain a result. 

Therefore, there is a need for faster way to determine the intersection of two 
polygons. 
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SUMMARY OF THE INVENTION 

To address these and other objectives, the present invention comprises a method 
for determining a polygonal intersection of a first polygon and a second polygon. An 
intersection of the boundary of the first polygon with the boundary of the second polygon 
is located by finding a node from which at least three portions (i.e., links) of polygon 
boundaries extend. From this node, a first portion of the boundary of the polygonal 
intersection is determined by identifying a portion of the boundary of the first polygon 
that is located inside the second polygon. Each subsequent portion of the boundary of the 
polygonal intersection is determined by selecting that portion of the boundary of either 
the first polygon or the second polygon that (1) connects to a leading end of a current 
portion of the boundary of the polygonal intersection and (2) forms the least angle with 
the current portion of the boundary of the polygonal intersection. 

According to a further aspect of the disclosed invention, intersected holes or 
islands are determined if they are in both polygons. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a diagram that illustrates a geographic region and data representations 
of some of the features located in the geographic region. 

Figure 2 is a diagram that illustrates a polygon used to represent a two- 
dimensional geographic area and a data representation of the polygon. 

Figure 3 is a block diagram illustrating a software program that determines the 
intersection of two polygons. 

Figure 4 depicts two intersecting polygons and illustrates operation of a process 
performed by the software program of Figure 3. 

Figures 5A-5C are a flow chart of the process performed by the software program 
of Figure 3. 

Figure 6 depicts a step in the process of Figures 5A-5C. 
Figures 7A and 7B depict another step in the process of Figures 5A-5C. 
Figure 8 depicts another step in the process of Figures 5A-5C. 
Figure 9 shows a polygon with a hole and data representation thereof. 
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1 Figure 10 shows an intersection of polygons with holes in the intersection. 

2 

3 DETAILED DESCRIPTION OF THE 

4 PRESENTLY PREFERRED EMBODIMENTS 

5 L THE GEOGRAPHIC DATABASE 

6 Referring to Figure 1, a geographic database 100 includes data that represent 

7 features that are located in a covered region 108. For example, the geographic database 

8 100 may contain data that represent roads 1 12 and points of interest 116. The geographic 

9 database 100 also includes data 120 that represent two dimensional features or entities 

10 124. Some of the different kinds of two dimensional features that can be represented in 

1 1 the geographic database 100 include land features (e.g., mountain ranges), bodies of 

12 water (e.g., lakes, ponds, swamps, wetlands), recreation areas (e.g., golf courses, parks, 

13 stadiums, forest preserves), business parks, industrial parks, and governmental 



€1 14 administrative areas (e.g., cities, towns, states, townships, municipalities, counties, school 



15 districts), as well as other types of areas, 
kjf 16 Referring to Figure 2, in the geographic database 100, two-dimensional features 

17 are represented using polygons. In other words, in the geographic database 100, the data 

18 that are used to represent two-dimensional features describe these features as polygons. 
P 19 In the geographic database 100, the boundary of the polygon corresponds to the boundary 



pi: 

s g; 20 of the represented two-dimensional feature. 



r? 21 The following terminology applies to the representation of two-dimensional 

22 features in the geographic database 100. 

23 "Node" - A point that terminates a link. 

24 "Line segment" - A straight line connecting two points. 

25 "Link" (or "edge") - A contiguous, non-branching string of one or more line 

26 segments terminating in a node at each end. 

27 "Shape point" - A point along a link between two nodes. 

28 "Oriented Link" - A link that has a starting node (referred to as the "reference 

29 node") and an ending node (referred to as the "non reference node"). 

30 "Simple Polygon" - An interior area of an outer boundary formed by a string of 

31 oriented links that begins and ends in one node. (A simple polygon does not cross itself.) 
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"Polygon" - The area bounded by an outer boundary and none or at least one 
interior boundary (e.g., a hole or island). A polygon is constructed from one outer simple 
polygon and none or at least one inner simple polygon. A polygon is simple if it just 
consists of one simple polygon, or complex if it has at least one inner simple polygon. 

The geographic database 100 follows certain conventions. For example, links do 
not cross themselves and not cross each other except at a node. Also, there are no 
duplicated shape points, nodes or links. Two links that connect to each other have a 
common node. In the geographic database 100, overlapping two-dimensional geographic 
features are represented by overlapping polygons. When polygons overlap, the boundary 
of one polygon crosses the boundary of the other polygon. In the geographic database 
100, the location at which the boundary of one polygon intersects the boundary of 
another polygon is represented by a node. (A node may be used to represent other 
locations along the boundary of a polygon other than a location at which the boundary of 
the polygon intersects the boundary of another polygon.) A shape point is not used to 
represent a point at which the boundary of a polygon intersects the boundary of another 
polygon. 

In the geographic database 100, a data representation 132 of a polygon that 
represents a two-dimensional feature includes a list 136 that identifies each oriented link 
that forms the boundary of the polygon. According to a convention used in the 
geographic database, the list 136 identifies the oriented links in clockwise order. If the 
polygon has a hole (i.e., a second polygon located entirely inside the first polygon), the 
hole is identified by a list of oriented links listed in counter clockwise order. 

II. METHOD FOR DETERMINING THE INTERSECTION OF POLYGONS 

As mentioned above, with applications that use geographic databases, there is 
sometimes a need to determine whether two two-dimensional features overlap and, if 
they do, to determine the two-dimensional area that forms the intersection (i.e., the 
overlapping portion) of the two two-dimensional areas. In terms of the polygons used to 
represent these two-dimensional features, there is a need to determine the intersection of 
polygons. 
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1 A method for determining the intersection of polygons is described in connection 

2 with Figures 3-10. Referring to Figure 3, the method is implemented by a software 

3 program 200 that uses the geographic database 100 that contains data representations of 

4 polygons, as described in connection with Figures 1 and 2. The software program 200 

5 receives two data representations of polygons as an input. (The two data representations 

6 of polygons conform to those described in Figure 2.) The software program 200 

7 provides, as an output, a data representation of each polygon that represents a polygonal 

8 intersection of the two data representations of polygons that were provided as input. 

9 Alternatively, the software program 200 provides an output indicating that there is no 

10 intersection of the two data representations of polygons that were provided as input. 

1 1 The software program 200 may be installed on a navigation system that uses a 
m 12 geographic database to provide navigation-related features to a driver of a vehicle. The 
%S 13 software program 200 may also be installed on a server that uses a geographic database to 
fjii 14 provide navigation-related services to various types of users, including vehicle drivers, 

^ 15 persons who are not drivers, users of personal digital assistants (PDAs), cell phone users, 

\H 16 and so on. The server that uses the geographic database may be connected to the Internet 

17 or other communications network. The software program 200 may also be used when 

18 compiling a geographic database, i.e., to produce a derived database product from a 

f=4: 19 source database. The software program 200 may also be used on a standalone computer 



IK 20 that uses a geographic database to provide various navigation-related, map-related, or 

H : 2 1 other types of features. 

22 In any of these different systems, the software program 200 may be stored on a 

23 computer-readable medium, loaded from the medium into a memory of the system and 

24 run, as needed. 

25 The determination of the polygonal intersection of two polygons is described with 

26 reference to Figures 4-10. Figure 4 shows two polygons, labeled P(A) and P(B). Using 

27 the disclosed embodiment, the intersection, P(I), of polygons P(A) and P(B) is 

28 determined. 

29 Figures 5 A-5C shows steps in a process 220 performed by the software program 

30 200 for determining the intersections of the polygons P(A) and P(B). The process 220 
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begins with a step in which the data representations of the two polygons are received 
(Step 224). The process 220 includes a step that forms minimum bounding rectangles for 
the two polygons P(A) and P(B) (Step 228). Figure 6 depicts formation of these 
minimum bounding rectangles for P(A) and P(B), labeled MBR(A) and MBR(B), 
respectively. The minimum bounding rectangles, MBR(A) and MBR(B), are compared 
to determine whether they intersect (Step 232). If the minimum bounding rectangles of 
P(A) and P(B) do not intersect, then the polygons P(A) and P(B) do not intersect, and a 
negative result is returned (Step 236). If the minimum bounding rectangles MBR(A) and 
MBR(B) intersect, then the polygons P(A) and P(B) may intersect and the process 220 
continues. 

The process continues by forming a node-link map. The node-link map is formed 
using the polygonal intersection, MBR(A B) of the minimum bounding rectangles 
MBR(A) and MBR(B) of the polygons P(A) and P(B). This polygonal intersection 
P(A_B) is determined as shown in Figure 6. To form the node-link map, each link of 
either polygon, P(A) or P(B), which is located entirely inside the polygonal intersection 
P(A_B) of the minimum bounding rectangles MBR(A) and MBR(B) of the polygons 
P(A) and P(B) is identified (Step 244). (A link is located entirely inside the intersection 
MBR(A_B) of the minimum bounding rectangles MBR(A) and MBR(B) if both ends of 
the link are inside the intersection MBR(A_B) of the minimum bounding rectangles 
MBR(A) and MBR(B).) Then, each node of each of these identified links is identified. 
(Shape points can be ignored at this step.) Each of these identified nodes is then included 
as an entry in the node-link map (Step 248). For each node included as an entry in the 
node-link map, all the links that connect to the node are identified and associated with the 
node in the map. 

The node-link map is used to identify a node (referred to herein as "N(START)") 
that is located on the polygon P(I) that forms the intersection of the input polygons P(A) 
and P(B). Using the node-link map, a node is identified that has three (or more) links 
associated with it (Step 252). A node from the node-link map that has three (or more) 
links associated with it is located on the polygon P(I) that forms the intersection of the 
two polygons P(A) and P(B). Because this node is located on the polygon P(I) that forms 
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1 the intersection of the two polygons P(A) and P(B), two (and only two) of the three or 

2 more links that connect to the node N(START) are links of the intersection polygon P(I). 

3 The next step is to identify which two of the three of more links form part of the 

4 intersection polygon P(I). 

5 In order to identify which of the links that connect to N(START) form part of the 

6 intersection polygon P(I), a process is performed that identifies one link of one of the 

7 original polygons, P(A) or P(B), that is located inside the other of the original polygons, 

8 P(B) or P(A) (Step 254). A link of one of the original polygons that is located inside the 

9 other of the original polygons is part of the intersection polygon P(I). This link is 

10 referred to as "LINK(START)." Details of the process for finding this link are illustrated 
in Figures 7 A, 7B and 8. 

Referring to Figures 7A and 7B Jn P(A) and P(B), all edges are oriented 

13 cj&ckwise. Node N is the starting node since four edges are connected to it, i.e., ael and 

14 ae2 from P(A), bel and be2 from ^(B). If ae2 is on the right side of bel and be2, then 
M 15 ae2 is inside P(B) and can be selected as the starting edge of the polygon intersection (as 
*II 16 shown in Figure 7A). If ael is/on the right side of bel and be2, then ael is inside P(B) 

17 and can be selected as the starting edge of the polygon intersection. If ae2 and ael both 

18 are not on the right side ofbel and be2, there is no intersection at N (as shown in 

H= 19 Figure 7B). To determine whether ae2 is on the right side of bel and be2, the process 

20 500 shown in Figure 8/can be used. Using the process 500 in Figure 8, the shape point 

21 bpl is selected fromybel, the shape point bp2 is selected from be2 and the shape point 

22 ap2 is selected from ae2. These shape points are selected because they are the closest 

23 shape points on tneir respective edges to node N(point p) although other shape points on 

24 these edges maybe selected as well. The process 500 in Figure 8 is used to test point if 

25 ap2 is on the right side of the line segment bpl_p and p_bp2. If so, ae2 is to the right of 

26 bel and be2 (as shown in Figure 7A). If not, ae2 is not on the right of bel and be2 (as 

27 shown in Figure 7B). 

28 After LINK(START) is identified, the next link that forms the intersection 

29 polygon P(I) is identified. This next link is identified by examining all the links that 

30 connect to the clockwise end of the LINK(START). Referring to Figure 5B, the 



-7- 



N0099US 





1 


clockwise end of the LINK(START) is determined from the clockwise direction of the 




2 


links in the original polygon (P(A) or P(B)) from which LINK(START) was taken 




3 


(Step 264). 




4 


If there is only one other link that connects to the clockwise end of the link, 




5 


LINK(START), then this link is the next link of the intersection polygon P(I). Data 




6 


indicating this next link are stored in the list formed to represent the intersection polygon 
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P(I) (Steps 272 and 276). 






tS/rf If there are two or more other links that connect to the clockwise end of the 




L^fNKCSTART), then the hrffthat forms the smallest angle in the counterclockwise 




10 


direction is the nexpifik of the intersection polygon P(I). Data indicating this next link 
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are stored in th€Tist formed to represent the intersection polygon P(I) (Steps 272, 278 and 
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276). S 
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After this link (i.e., the link after LINK(START)) is added to the list formed to 
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represent the intersection polygon P(I), the process continues by determining the next 
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link that forms the intersection polygon P(I). This next link is determined in a similar 




16 


manner as just described. That is, all the links that connect to the clockwise end of the 




17 


link are identified (Steps 264 and 272, again). If there is only one other link, then this 


u 
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18 


link is added to the list of links that form the intersection polygon P(I). If there are two or 




19 


more other links, then the next link is the one that forms the smallest angle in a 
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20 


counterclockwise direction from the link most recently identified as being part of the 




21 


intersection polygon P(I) (Step 278, again). 




22 


The process continues until the N(START) is encountered (Step 280). When 




23 


N(START) is encountered, then one simple polygon (that forms part or all of the entire 




24 


polygon intersection P(I) of P(A) and P(B)) has been identified. Data describing this 




25 


simple polygon are contained in the list 258. These data are added to a simple polygon 




26 


list in a data representation 318 (in Figure 5C) of the entire intersection polygon P(I) of 




27 


P(A) and P(B) (Step 284). 




28 


All used links and nodes in the link-node map are removed (Step 288). If there 




29 


are still nodes with three or more connected links in the node-link map after a simple 




30 


polygonal intersection has been identified (Step 292 in Figure 5B), then there is at least 
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1 one more distinct simple polygonal intersection between the first polygon and the second 

2 polygon. The process is repeated (starting at Step 252 in Figure 5A) to determine each 

3 additional simple polygonal intersection. The process is repeated until there is no node in 

4 the link-node map that has more than three links connected to it. 

5 The process 220 finds intersected holes and islands in P(A) and P(B). Some two- 

6 dimensional features represented in the geographic database 100 are encompassed 

7 entirely within other represented two-dimensional features. As an example, an island is 

8 located entirely within a body of water. In the geographic database 100, a two- 

9 dimensional feature that is encompassed entirely within another two-dimensional is 

10 represented as a hole inside a polygon. 

1 1 The data representation of the polygon containing a hole includes (1) a list of 

p~ 12 edges that define the boundary of the entire polygon (i.e., similar to Figure 2, described 

O 13 above) and (2) a list of edges that define the boundary of the hole located inside the 

14 polygon. According to a convention used in the geographic database 100, the list of 

15 edges in the data representation of a hole inside a polygon are ordered in counter- 

ijfli 16 clockwise order. Figure 9 shows a polygon with a hole and data representation thereof. 

? 17 When determining the intersection of two polygons that have holes, it is possible 

y 18 that some of the holes are included in the intersection. Figure 10 illustrates the 

h& 19 intersection of two polygons having holes/islands in their intersection. 

20 When the disclosed process 220 is used to determine the polygonal intersection of 

21 two polygons, the data representation of the polygonal intersection includes a 

22 representation of any holes of either of the original polygons that are located inside the 

23 polygonal intersection. First, a determination is made whether the minimum bounding 

24 rectangle of the hole is located inside the intersection of the minimum bounding 

25 rectangles of the two original polygons. If the minimum bounding rectangle of the hole 

26 is not located completely inside the intersection of the minimum bounding rectangles of 

27 the two original polygons, then the hole is not located in the polygonal intersection. If 

28 the minimum bounding rectangle of the hole is located inside the intersection of the 

29 minimum bounding rectangles of the two original polygons, then a determination is made 

30 whether any point (e.g., any node) of the hole is located inside the other of the polygons. 
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1 Figure 5C shows steps in the process 220 for handling holes/islands in the 

2 intersection polygon, P(I). Starting with an island/hole from polygon A, a search is made 

3 to find an immediate parent simple polygon in polygon B. If there is a simple polygon in 

4 B and its orientation is clockwise, then the island/hole from polygon A is a simple 

5 intersection polygon. If there is no parent polygon in B or the orientation of the parent 

6 simple polygon in B is counter clockwise, then the island/hole from polygon A is not a 

7 simple intersection polygon (Step 302). These steps are repeated until all islands/holes 

8 from polygon A have been considered. 

9 The process continues by considering an island/hole from polygon B and 

10 searching for an immediate parent simple polygon in polygon A. If there is a simple 

1 1 polygon in A and its orientation is clockwise, then the island/hole from polygon B is a 

12 simple intersection polygon. If there is no parent polygon in A or the orientation of the 

13 parent simple polygon in A is counter clockwise, then the island/hole from polygon B is 
fil 14 not a simple intersection polygon (Step 306). This process is repeated until all 

I . § 15 islands/holes from polygon B have been visited. 

^ 16 After accounting for any holes or islands in the intersection, a data representation 

M 

ii 17 of the entire complex polygon intersection, P(I), is built from the list(s) of simple 

□ 

5t 18 polygons list 318 (Step 307). The intersection of P(A) and P(B) may include one or more 

K- 19 polygons, each of which may be either simple or complex. The data representation of the 

P 20 polygonal intersection can then used by other applications in the system in which the 

21 software program 200 is installed that require data indicating the polygonal intersection 

22 of the P(A) and P(B) (Step 288) 
23 

24 in. ADVANTAGES 

25 The present method for determining the intersection of polygons is 

26 computationally less intensive than prior methods. The present method achieves this 

27 advantage, in part, because only the polygon edge is considered instead of entire line 

28 segments, therefore greatly reducing the amount of computation. 
29 
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1 IV. FURTHER CONSIDERATIONS 

2 Consistent ordering - clockwise or counterclockwise 

3 In the embodiments described above, it was explained that the data 

4 representations of polygons included lists of links used to represent the boundaries of the 

5 polygons and that the links contained on a list were ordered corresponding to a clockwise 

6 traversal of the represented polygon boundary via the links. It is not required that the 

7 lists of links used to represent polygons identify the links in clockwise order so long as 

8 the order in which the links are listed is consistent throughout the geographic database. 

9 For example, according to an alternative embodiment, instead of listing the links that 

10 form polygon boundaries in a clockwise order, data representations of polygons could list 

1 1 the links in a counterclockwise order. If a counterclockwise order is used for polygon 

12 boundaries, it might be preferable to use a clockwise order to represent holes in polygons. 
S> 13 

flj 14 It is intended that the foregoing detailed description be regarded as illustrative 

f 7* : 15 rather than limiting and that it is understood that the following claims including all 

& 16 equivalents are intended to define the scope of the invention. 
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